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Abstract 

We have searched for the lepton flavor violating decays r~ — > £~K§ (t = e or /x), using a data 
sample of 281 fb _1 collected with the Belle detector at the KEKB e + e~ asymmetric-energy collider. 
No evidence for a signal was found in either of the decay modes, and we set the following upper 
limits for the branching fractions: B(t~ — > e~K®) < 5.6 x 10~ 8 and B(t~ — > u~K§) < 4.9 x 10 -8 
at the 90% confidence level. These results improve the previously published limits set by the CLEO 
collaboration by factors of 16 and 19, respectively. 

PACS numbers: 11.30.Fs; 13.35.Dx; 14.60.Fg 
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INTRODUCTION 



Lepton flavor violation (LFV) is allowed in many extensions of the Standard Model (SM), 
such as Supersymmetry (SUSY) and leptoquark models. In particular, lepton flavor violating 
decays with K$ mesons are discussed in models with heavy singlet Dirac neutrinos 0, 
R— parity violation in SUSY 0,111 , dimension-six effective fermionic operators that induce 
t — n mixing p[. Experiments at the I?-factories allow searches for lepton flavor violating 
decays with a very high sensitivity. The best upper limits of Bir~ — > e~K®) < 9.1 x 10~ 7 and 
B(t~ — > fi~Kg) < 9.5 x 10~ 7 at the 90% confidence level were set by the CLEO experiment 
using 13.9 fb _1 of data (if. 

In this paper, we report a search for the lepton flavor violating decays r~ — * £~K§ 
(£ = e or /i)[f] using 281 fb _1 of data collected at the T(4S) resonance and 60 MeV below 
it with the Belle detector at the KEKB e + e~ asymmetric-energy collider jfj. 

The Belle detector is a large-solid-angle magnetic spectrometer that consists of a silicon 
vertex detector (SVD), a 50-layer central drift chamber (CDC), an array of aerogel threshold 
Cerenkov counters (ACC), a barrel-like arrangement of time-of-flight scintillation counters 
(TOF), and an electromagnetic calorimeter comprised of CsI(Tl) crystals (ECL), all located 
inside a superconducting solenoid coil that provides a 1.5 T magnetic field. An iron flux- 
return located outside of the coil is instrumented to detect K® mesons and to identify muons 
(KLM). The detector is described in detail elsewhere 0]. 

Particle identification is very important in this measurement. We use particle identi- 
fication likelihood variables based on the ratio of the energy deposited in the ECL to the 
momentum measured in the SVD and CDC, the shower shape in the ECL, the particle range 
in the KLM, the hit information from the ACC, the measured dE/dX in the CDC and the 
particle's time-of-flight from the TOF. For lepton identification, we form a likelihood ratio 
based on the electron probability V(e) [|[ and the muon probability V(\i) determined by 
the responses of the appropriate subdetectors. 

For Monte Carlo (MC) simulation studies, the following programs have been used to 



generate background events: KORALB/TAUOLA 11 OB for r + r , QQ |ll| for BB and con- 
tinuum, BHLUMI H3 for Bhabha events, KKMC [l| for e+e" -> fj+fj,- and AAFH 
for two-photon processes. Since the QQ generator does not include some rare processes 
that potentially contribute to final states with a K® meson, we generated special samples of 
e + e~ — > D* + D^~, a process that was recently observed by the Belle collaboration [15j- Sig- 
nal MC is generated by KORALB/TAUOLA. Signal r decays are two-body and assumed to 
have a uniform angular distribution in the r lepton's rest frame. The Belle detector response 
is simulated by a GEANT 3 0] based program. All kinematic variables are calculated in 
the laboratory frame unless otherwise specified. In particular, variables calculated in the 
e + e~ center-of-mass (CM) frame are indicated by the superscript "CM" . 

DATA ANALYSIS 

We search for t + t~ events in which one r (signal side) decays into £K§ (Kg — > 7r + 7r~), 
while the other r (tag side) decays into one charged track (with a sign opposite to that 
of the signal-side lepton) and any number of additional photons and neutrinos. Thus, the 



[|] Unless otherwise stated, charge conjugate decays are included throughout this paper. 
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experimental signature is: 



{t~ -> £-(= e" or fi~) + Kg(^> tt+tt")} + {r+ -> (a track) + + (n TAG 7 > 0) + A" (missing)}. 

All charged tracks and photons are required to be reconstructed within a fiducial volume, 
defined by —0.866 < cos# < 0.956, where 9 is the polar angle with respect to the direction 
opposite to the e + beam. We select charged tracks with momenta transverse to the e + beam 
p t > 0.1 GeV/c and photons with energies E 1 > 0.1 GeV. 

Candidate r-pair events are required to have four charged tracks with a zero net charge. 
Events are separated into two hemispheres corresponding to the signal (three-prong) and 
tag (one-prong) sides by the plane perpendicular to the thrust axis The magnitude 

of the thrust is required to be larger than 0.9 to suppress the qq continuum background. 
The Kg is reconstructed from two oppositely-charged tracks in the signal side that have an 
invariant mass 0.482 GeV/c 2 < M n + n - < 0.514 GeV/c 2 , assuming a pion mass for both 
tracks. The 7r + 7r~ vertex is required to be displaced from the interaction point (IP) in the 
direction of the pion pair momentum ^| . In order to avoid fake Kg candidates from photon 
conversions (i.e. 7 — > e + e~), the invariant mass reconstructed by assigning the electron mass 
to the tracks, is required to be greater than 0.2 GeV/c 2 . The signal side track not used in 
the K$ reconstruction is required to satisfy the lepton identification selection. The electron 
and muon identification criteria are V(e) > 0.9 with p > 0.3 GeV/c and V{ji) > 0.9 with 
p > 0.6 GeV/c, respectively. After the event selection described above, most of the remaining 
background comes from generic r + r~ and continuum events that contain a real Kg meson. 

To ensure that the missing particles are neutrinos rather than photons or charged particles 
that fall outside the detector acceptance, we impose additional requirements on the missing 
momentum vector, p m i SS , calculated by subtracting the vector sum of the momenta of all 
tracks and photons from the sum of the e + and e~ beam momenta. We require that the 
magnitude of p m i SS be greater than 0.4 GeV/c and that its direction point into the fiducial 
volume of the detector, as shown for the — ► fi~Kg mode in Fig. ^ (a) and (b). The 
total visible energy in the CM frame, E™, is defined as the sum of the energies of the 
Kg candidate, the lepton, the tag-side track (with a pion mass hypothesis) and all photon 
candidates. We require E™ to satisfy the condition: 5.29 GeV < E™ < 10.0 GeV (see 
Fig. [U(c)). Since neutrinos are emitted only on the tag side, the direction of p m i SS should 
lie within the tag side of the event. The cosine of the opening angle between p m i SS and 
the tag-side track in the CM system, cos0^_ miss , is therefore required to be greater than 
(see Fig. Q (d))- For all kinematic distributions shown in Fig. reasonable agreement 
between the data and background MC is observed. In order to suppress background from 
qq (q — u,d,s, c) continuum events, the following requirements on the number of the photon 
candidates on the signal and tag side are imposed: n SIG < 1 and n TAG < 2, respectively. 

Finally, the correlation between the reconstructed momentum of the £Kg system, ptx s , 
and the cosine of the opening angle between the lepton and Kg, cos0£k s is employed to 
further suppress background from generic t + t~ and continuum events via the requirements: 
cos6tK s < 0.14 x \og{p iKs - 2.7) + 0.7, where p iKs is in GeV/c (see Fig. EJ). While this 
condition retains 99% of the signal, 99% of the generic t + t~ and 84% of the uds continuum 
background are removed. Following all the selection criteria, the signal detection efficiencies 
for the t~ — > e~Kg and t" — > [i~Kg modes are 15.0% and 16.2%, respectively. 
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FIG. 1: Kinematic distributions used in the event selection after K§ mass and muon identification 
requirements: (a) the momentum of the missing particle; (b) the polar angle of the missing particle; 
(c) the total visible energy in the CM frame; (d) the opening angle between the missing particle and 
tag-side track in the CM frame. The signal MC distributions are indicated by the filled histograms, 
the total background including t + t~ and qq is shown by the open histogram, and solid circles are 
data. While the signal MC (r _ — > fi~K§) distribution is normalized arbitrarily, the data and 
background MC are normalized to the same luminosity. Selected regions are indicated by arrows 
from the marked cut boundaries. 



RESULTS 

Signal candidates are examined in the two-dimensional plots of the £~K§ invariant mass, 
M £K o, and the difference of their energy from the beam energy in the CM system, AE. A 
signal event should have M iK o close to the r-lepton mass and AE close to 0. For both modes, 
the M eK o and AE resolutions are parameterized from the MC distributions around the peak 
with bifurcated Gaussian shapes to account for initial state radiation. These Gaussian have 
widths a^ gh/ low = 6.2/7.4 MeV/c 2 and a^| h/low = 20/26 MeV for the r" -> e~K% mode, 

and a^ gh/ low = 6.1/5.9 MeV/c 2 and a^| h/low = 19/23 MeV for the r~ -> yTK% mode, where 

the "high/low" superscript indicates the higher/lower side of the peak. 

We blind a region of ±5om around the r mass in M£ K o and a region of —0.5 GeV < 

eK s s 
AE < 0.5 GeV so as not to bias our choice of selection criteria. Figure |3] shows scatter-plots 

for data and signal MC samples distributed over ±15<r in the M eK o — AE plane. Most of 
the surviving background events in both modes come from — > -k^K^ and D — > i^vK^ 
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FIG. 2: Scatter-plots of (a) signal MC (r" -> n'Kg), (b) continuum MC, (c) generic r+r" MC 
events and (d) data on the p^K s vs cos 0^k s plane. Selected regions are indicated by curves with 
arrows. 



decays. The remaining continuum backgrounds in the r~ — > ^~K^ mode are combinations 
of a true Kg meson and a fake lepton. 

To optimize our search sensitivity, we select an elliptically shaped signal region of mini- 
mum area with the same signal acceptance as that of a rectangular box corresponding to ±5<r 
in the MC resolution for the M eK o — AE plane. The signal efficiencies after all requirements 
are 11.8% for the r~ — > e~K® and 13.5% for the r~ — > n~K$, respectively. 

As there are few remaining MC background events in the signal ellipse, we estimate 
the background contribution using the M eK o sideband regions defined by rectangular areas 
beside the signal ellipse shown in Fig. |H] (a) and (b). Extrapolation to the signal region 
assumes that the background distribution is flat in M tK a. We find the expected background 
in the ellipse to be 0.2 ± 0.2 events for both modes. Finally, we uncover the blinded region 
and find no data events in the signal region of the r~ — > e~K® and r" — > f^~Kg modes 
(see Fig. El (a) and (b)). Since no statistically significant excess of data over the expected 
background in the signal region is observed, we apply a frequentist approach to calculate 
upper limits on the signal yields 0|. The resulting limits for the signal yields at 90% 
confidence level, S90, are 2.23 events in both modes. The upper limits on the branching 
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FIG. 3: Scatter-plots of data in the M eK o - AE plane: (a) and (b) correspond to the ±15cr 
area for the r~ — ► e~Kg and r~ — ► fJ-^K^ modes, respectively. The data are indicated by the 
solid circles. The filled boxes show the MC signal distribution with arbitrary normalization. The 
elliptical signal region shown by a solid curve is used for evaluating the signal yield. The vertical 
dashed lines denote the boundaries of the blind regions, while the rectangular areas beside the 
signal region are used to estimate the expected background in the elliptical region. 



B(r~ i-Ki) < o ^ O90 , + ,_ w (1) 



fraction before the inclusion of systematic uncertainties are then calculated as 

2eB(K£ ^+ix~)N T 

where B(Kl -> tt+tt-) = 0.6895 ± 0.0014 and iV TT = 251 x 10 6 is the number of r-pairs 
produced in 281 fb _1 of data. We obtain N TT using a TT = 0.892±0.002 nb, the e + e~ — > t + t~ 
cross section at the T(4S) resonance calculated by KKMC 0]. The resulting values are 
B{t~ -> e-Kg) < 5.5 x 10~ 8 and B{r~ -> (irKg) < 4.8 x 10" 8 . 

The dominant systematic uncertainties on the detection sensitivity: 2eN TT B(K® — > 
7r + 7r~) come from K$ reconstruction and tracking efficiencies. These are 4.5% and 4.0%, 
respectively, for both modes. Other sources of the systematic uncertainties are: the trigger 
efficiency (0.5%), lepton identification (2.0%), MC statistics (0.3%), branching fraction of 
Kg — > 7r + 7r~ (0.2%) and luminosity (1.4%). Assuming no correlation between them, all 
these uncertainties are combined in quadrature to give a total of 6.5%. 

While the angular distribution of r~ — * £~K§ decay is initially assumed to be uniform in 



this analysis, it is sensitive to the lepton flavor violating interaction structure [21(. The spin 



correlation between the r lepton in the signal and that in the tag side must be considered. 
A possible nonuniformity was taken into account by comparing the uniform case with those 
assuming V — A and V + A interactions, which result in the maximum possible variations. 
No statistically significant difference in the M eK o - AE distribution or the efficiencies is 
found compared to the case of the uniform distribution. Therefore, systematic uncertainties 
due to these effects are neglected in the upper limit evaluation. 
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Upper limits on the branching fractions at the 90% confidence level including these sys- 
tematic uncertainties are calculated with the POLE program without conditioning 2^ . 
The resulting upper limits on the branching fractions at the 90% confidence level are 

B(t~ -> e'Kl) < 5.6 x 10~ 8 
B(r~ -> p-Kg) < 4.9 x 10~ 8 . 



DISCUSSION 

In the R— parity violating SUSY scenario, there are three kinds of terms (A, A' and A") 
with a total of 45 couplings. In this model, r~ could decay into l~ K$ via tree- level scalar 
neutrino exchange by the AA' couplings. Using our results, the limits on the products AA' 
as a function of the scalar neutrino mass (Ma) are given as 



\X i31 X' a2 \(i = 1,2), \X i31 X' m \(i = 2,3) < 4.5 x 10- 4 (il4/100GeV/c 2 ) 2 for r" -> e"^ 



|A i32 A' il2 |(z = 1,2), |A i23 A; 21 |(2 = 1,3), < 4.1 x 10- 4 (M,>/100GeV/ C 2 ) 2 for r" -> fi~ K° 



where i is the generation number. These bounds are more stringent than the previous bounds 
obtained in R— parity violating models from r~ decay including a pseudoscalar meson 0, HI . 

The improved sensitivity to rare r lepton decays achieved in this work can be used to 
constrain the new physics scale for the dimension-six fermionic effective operators involving 
t — fi flavor violation, motivated by neutrino oscillations 0]. From our upper limit for the 
branching fraction of the r~ —* n~~K§ decay, lower bounds of 36.2 TeV and 37.2 TeV can 
be obtained for the axial-vector and pseudoscalar operators, respectively. 



CONCLUSION 

In conclusion, we have searched for the lepton flavor violating decays r~ — > £~K§ (£ = 
e or fi) using data collected with the Belle detector at the KEKB e + e~ asymmetric-energy 
collider. We found no signal in either mode. The following upper limits on the branching 
fractions at the 90% confidence level are obtained: B(r~ — > e~K§) < 5.6 x 10~ 8 and 
B(r~ — > fi~Kg) < 4.9 x 10~ 8 . These results improve the search sensitivity by factors of 16 
and 19, respectively, compared to the previous limits obtained by the CLEO experiment. 



Acknowledgments 

We thank the KEKB group for the excellent operation of the accelerator, the KEK cryo- 
genics group for the efficient operation of the solenoid, and the KEK computer group and 
the National Institute of Informatics for valuable computing and Super-SINET network sup- 
port. We acknowledge support from the Ministry of Education, Culture, Sports, Science, 
and Technology of Japan and the Japan Society for the Promotion of Science; the Australian 
Research Council and the Australian Department of Education, Science and Training; the 
National Science Foundation of China and the Knowledge Innovation Program of Chinese 
Academy of Sciencies under contract No. 10575109 and IHEP-U-503; the Department of Sci- 
ence and Technology of India; the BK21 program of the Ministry of Education of Korea, and 
the CHEP SRC program and Basic Research program (grant No. R01-2005-000-10089-0) of 



8 



the Korea Science and Engineering Foundation; the Polish State Committee for Scientific 
Research under contract No. 2P03B 01324; the Ministry of Science and Technology of the 
Russian Federation; the Slovenian Research Agency; the Swiss National Science Founda- 
tion; the National Science Council and the Ministry of Education of Taiwan; and the U.S. 
Department of Energy. 



[1] A. Ilakovac, Phys. Rev. D 62, 036010 (2000). 

[2] J. P. Saha and A. Kundu, Phys. Rev. D 66, 054021 (2002). 

[3] R. Barbier et al, Phys. Rep. 420, 1 (2005). 

[4] D. Black et al, Phys. Rev. D 66, 053002 (2002). 

[5] S. Chen et al. (CLEO Collaboration), Phys. Rev. D 66, 071101 (2002). 

[6] S. Kurokawa and E. Kikutani, Nucl. Instr. and Meth. A 499, 1 (2003), and other papers 

included in this Volume. 

[7] A. Abashian et al. (Belle Collaboration), Nucl. Instr. and Meth. A 479, 117 (2002). 

[8] K. Hanagaki et al., Nucl. Instr. and Meth. A 485, 490 (2002). 

[9] A. Abashian et al., Nucl. Instr. and Meth. A 491, 69 (2002). 

[10] S. Jadach and Z. Was, Comp. Phys. Commun. 85, 453 (1995). 

[11] QQ is an event generator developed by the CLEO Collaboration and described in 
http://www.lns.cornell.edu/public/CLEO/soft/QQ/. It is based on the LUND Monte Carlo for 
jet fragmentation and e + e~ physics described in T. Sjostrand, Comp. Phys. Commun. 39, 347 
(1986) and T. SjSstrand, Comp. Phys. Commun. 43, 367 (1987). 

[12] S. Jadach et al, Comp. Phys. Commun. 79, 305 (1992). 

[13] S. Jadach et al, Comp. Phys. Commun. 130, 260 (2000). 

[14] F. A. Berends et al, Comp. Phys. Commun. 40, 285 (1986). 

[15] T. Uglov et al (Belle Collaboration), Phys. Rev. D 70, 071101 (2004). 

[16] R. Brun et al, GEANT 3.21 CERN Report No. DD/EE/84-1, 453. 

[17] S. Brandt et al, Phys. Lett. 12, 57 (1964); E. Farhi, Phys. Rev. Lett. 39, 1587 (1977). 

[18] K. Sumisawa et al. (Belle Collaboration), Phys. Rev. Lett. 95, 061801 (2005). 

[19] G.J. Feldman and R.D. Cousins, Phys. Rev. D 57, 3873 (1998). 

[20] S. Eidelman et al. (Particle Data Group), Phys. Lett. B 592, 1 (2004). 

[21] R. Kitano and T. Okada, Phys. Rev. D 63, 3873 (2001). 

[22] J. Conrad et al, Phys. Rev. D 67, 012002 (2003). 



9 



